Far infrared cyclotron resonance and Faraday effect in Bi2Se3 
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The complex Faraday angle and transmission for the topological insulator Bi2Se3 were measured 
in the far infrared in magnetic fields up to 8 T and analyzed in terms of a simple Drude-Lorentz 
magneto-plasma dielectric model. Bulk carriers dominate the optical response. The bulk carriers 
are electrons as we determined from the sign of the Faraday angle. We obtain the bulk band edge 
cyclotron mass mcb = 0.16 ± O.Olme and free electron concentration in the 10^^ range. Electron- 
phonon interaction effects were found to be weak. 
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PACS numbers: 78.20.Ls 03.65.Vf 

I. INTRODUCTION 

It was shown recently theoretically^'^ and confirmed 
experimentally by ARPES^ii that Bi2Se3 and Bi2Te3 
have an odd number (one) of Dirac cone surface states in 
the unit cell and thus are nontrivial topological insula- 
tors. This makes them distinct from graphene which also 
has Dirac electrons but an even (two) number of cones 
per unit cell. Predictions of several novel properties have 
generated a great interest in these materials ^i"— 

Although the optical properties of Bi2Se3 have been 
studied in the past, only recently it has been probed for 
unusual far infrared magneto-optical response.™ We have 
grown Bi2Se3 crystals with low free electron concentra- 
tion which enabled observation of the cyclotron resonance 
and Faraday rotation below the strong low frequency po- 
lar phonon. In this paper, we report results of these 
measurements together with the model calculations. 



II. EXPERIMENTAL 

Single crystals of Bi2Se3 were grown as described 
in Ref. @, and we used a sample from batch iv. A 
2x2x0.014 mm flake was exfoliated and flattened on a 
flrm surface. The sample thickness was obtained from 
etalon oscillations and confirmed by electron microscope 
images of the sliced flake. Infrared data were obtained 
on a Bomem DAS Fourier-transform (FTIR) spectrome- 
ter and an optically pumped far infrared laser at a fre- 
quency of 42.3 cm~^. The samples were placed either 
in an optical cryostat for B = measurements or in an 
optical split-coil magnet for measurements in magnetic 
fields up to 8 T in the Faraday geometry: kughtWB and 
-B||c-axis. Our technique for the complex Faraday angle 
measurements in the far infrared is described in Refs. 
andli. 



III. MODEL 

The experiments reported here differ from conventional 
cyclotron resonance in semiconductors because the rela- 



tively high carrier density and low lying phonon leads 
to strong magneto-plasma effects. This results in a very 
low transmission near the plasma frequency and multiple 
reflections inside the flake in transparent regions (etalon 
effect) which produces a complicated optical response. 
Therefore, it is necessary to model the magneto-optical 
response in order to obtain the material parameters. We 
model our experiment as transmission and reflection of a 
slab with bulk optical constants and (possibly) two iden- 
tical 2D layers with surface optical constants on both sur- 
faces of the slab. For our free carrier concentration level, 
the Fermi energy is less than 20 meV above the bottom of 
the lowest conduction band. This is the only filled band 
because, according to Ref.'o', the second conduction band 
minimum is 40 meV higher. Since huj and Ep are small 
compared with the band gap ~300 meV, a parabolic dis- 
persion of the bulk electrons is assumed. Therefore, ex- 
periment gives essentially the band edge mass. 

In this case, we can use a Drude magneto-optical re- 
sponse of the free electron gas (1st term in Eq. [T]). The 
corresponding complex dielectric function for the bulk is: 
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where j enumerates the phonon oscillators, ujj is the TO 
phonon frequency, Sj is the phonon spectral weight, a 
subscript b refers to bulk electrons, is the dielec- 
tric constant at higher infrared frequencies, ujpb is the 
unscreened plasma frequency for free electrons: w^j, = 
An N^e^ / ml , where Ni,, e and are concentration, 
charge, and an effective mass of electrons; 7b and are 
the damping rates; uJcb — eB/{mcbc) is the cyclotron 
frequency, where B is magnetic field, met is bulk cy- 
clotron mass, and c is speed of light. The ± sign refers 
to the right and left circular polarized electromagnetic 
modes. For a fixed (laser) frequency, we use a quarter- 
wave plate to obtain circular polarized light. In case of 
the far infrared spectra, we used linear polarized light 
thus measuring an average of two circular polarizations. 
For a general polarization state, measured transmission 
is T = (0.5-l-;3)r+-|-(0.5-/3)T-, where -0.5 < /3 < 0.5. 
The surface states of a topological insulator are 2D 
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Dirac quasiparticles that split into Landau levels in an 
applied B\\c magnetic field : 
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where integer n enumerates energy levels, e and vp are 
charge and Fermi velocity of the Dirac particles. The full 
band dispersion for Bi2Se3 has been recently reported in 
Ref. llOi The optical dipole selection rules allow tran- 
sitions |n| — |n| -|- 1 for the cyclotron resonance active 
mode of circularly polarized light. The linear dispersion 
leads to the square root dependence on n and B. How- 
ever, in our sample the chemical potential /i is ~ 10 meV 
above the bottom of the conduction band which, in turn, 
is ~ 200 mcV above the Dirac point. Therefore, we find 
n w 25 for 8 T magnetic field. At these values of n, we 
can use an approximate formula 
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and the cyclotron resonance becomes semiclassical. The 
combination En/vp plays the role of the cyclotron mass 
rrics for surface states. Therefore, we model the optical 
response of the surface states as a 2D free electron gas. 

We use conventional expressions for transmission and 
refiection of a slab (medium 2, bulk) between media 1 
and 3 (both vacuum here) for light propagating from 1 
to 3 with the following modifications. The presence of 
surface states at the interfaces 1-2 and 2-3 changes the 
complex Fresnel transmission and reflection coefficients: 
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where ni = — 1, n2 — ^Je^Juj) are refractive indices, 
and surface admittance is 



Wo 



(5) 



where ujps = 47rA^s e^/(m*c) is a 2D plasma frequency, a 
subscript s refers to the surface states. Using standard 
expressions for the total complex transmission t and re- 
flection r amplitude coefficients of a slab 
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where (j) — exp (i27m2 Lod) and d is a sample thickness, we 
calculate transmission = \t^\'^n^/ni and refiection 
= \r^\^ intensity coefficients. We define the complex 
Faraday angle Op via tan^i? = tyx/txx, where standard 



TABLE I. Fit parameters for Bi2Se3 flake and calculated val- 
ues (last 3 columns). 



Set" 
(#" 


Wp6 

(cm-i) 


lb 
(cm-i) 


(me) 


iocb (8 T) 
(cm-i) 


(10^^ cm^^*) 


(meA 


1 (FTIR) 


382 


8 


0.15 


50 


2.3 


8 


2-1 (laser) 


385 


7.5 


0.16 


48 


2.6 


8 


2-2 (laser) 


750 


7.5 


0.16 


48 


9.9 


20 


3 (FTIR) 


465 


8 


0.16 


48 


3.8 


11 


4 (laser) 


520 


8 


0.17 


44 


5.1 


12 



In chronological order. Fixed phonon parameters are 
u)i = 67 cm-l, Si = 3.94 X 10^ cm-^, 71 = 5 cm-l, 
ij2 = 134 cm-i, 52 = 7056 cm-^, 72 = 2 cm"^, and Co 



: 25.6. 



m. 



transformation between circular and linear polarization 
bases = txx ± itxy is assumed. 

We estimate Ups using the ARPES data^*^ and ex- 
pressions for 2D free Fermi gas (one spin per state): 
For Ep K, 10 meV, the surface states are 



below kp = 0.08 A^^ which gives Ng 

^ assuming to* — O.ISto, 



5 X 10 cm 



and ujps ~ 20 cm 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 



Zero field far infrared results for this sample were re- 
ported earlier Set #1 shown in Table 1 is the result 
of a simultaneous fit of zero field reflectivity and trans- 
mission. Other sets of fit parameters shown in Table 1 
model results from different magnetic field experiments 
discussed below. Data of experiment #2 were fitted by 
the model of two crystals with parameters #2-1 and #2- 
2 and effective areas 82 % and 18 %, correspondingly. 
We understand variation of the bulk concentrations from 
one experiment to another and within experiment #2 
as a consequence of nonuniform carrier concentration of 
the sample. Characteristic size of the regions with differ- 
ent concentrations should be on the order of the infrared 
laser focal spot (« 1 mm). In FTIR measurements, the 
light spot was larger than the sample aperture. We also 
cannot rule out a possible increase of the effective bulk 
concentration with time due, for example, to exposure of 
the flake to air during the few months separating mea- 
surements #1 and #4. 

Figure [T] shows the measured and calculated trans- 
mission spectra of the Bi2Se3 flake. Oscillations of the 
transmission in zero field appearing below 60 cm~^ are 
due to the etalon effect within the transparency band 
between the plasma frequency and the TO phonon fre- 
quency. In non zero fields, the magneto-plasma dielec- 
tric function (1) changes due to the cyclotron resonance, 
and the etalon features for left and right circular polar- 
ization split in frequency and change their amplitudes. 
The material becomes more transparent in magnetic field 
because the T"*" (cyclotron resonance inactive) polariza- 
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FIG. 1. (Color online) Transmission of Bi2Se3 flake for un- 
polarized light in magnetic field B\\c: (a) measured and (b) 
calculated; green dash-dot curve is calculated with surface 
states (SS) included. 



tion becomes less absorbing. In the fit ^^3 shown in 
fig. [ijb) , we match frequency positions of the main max- 
ima in the experimental and model transmissions. The 
green dash-dot curve in fig. Wijo) shows how much the 
transmission changes if surface states are added with an 
exaggerated strength uips = 100 cm~^ (7^ = 8 cm~^, 
rrics = iTicb = 0.16TOe). The curve calculated for the 
estimated ARPES value of ujps = 20 cm"^ is almost in- 
distinguishable from the ujps — red dashed curve. 

We have also measured unpolarized magneto- 
transmission of the sample at the fixed laser frequency 
42.3 cm""'^. Figure [2] shows measured and calculated 
transmission spectra of the same Bi2Se3 sample as a 
function of applied magnetic field. The upturn of the 
experimental curve above 5 T cannot be reproduced by 
our simple model with one concentration either with or 
without surface states. However, we can fit the entire 
experimental curve by the weighted average of two 
independent sets of bulk parameters, simulating inho- 
mogeneity of the free carrier concentration, described 
by sets #2-1 and 2-2. Including different thicknesses 
associated with each set of bulk parameters to concur- 
rently simulate thickness variations can not significantly 
reduce the concentration difference between sets #2-1 
and 2-2 if the thicknesses are limited to 14 ± 4 microns 
as measured by scanning electron microscope. Based 
upon these simple models, it appears that concentration 
inhomogeneity may cause deviations from the simple 
line shape expected from a single uniform bulk response. 
A distribution of concentrations will cause a distribution 
of plasmons (zero crossings of the real part of epsilon) 
which will affect the optical response in a nontrivial way. 
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FIG. 2. (Color online) Measured (blue solid lines) and cal- 
culated (bulk only) transmission of nearly unpolarized light 
(fi = —0.02) for Bi2Se3 flake at the flxed laser frequency 
42.3 cm~^. (a) The calculated transmission for the regions of 
the flake with low carrier concentration (green dotted curve), 
high concentration (red dash-dotted curve), and a weighted 
admixture of the two which flts the 10 K data (pink dashed 
curve), (b) The same model and parameters used in (a) 
flt the 50 K data provided the scattering rate is raised to 
71 = 72 = 13.5 cm~^ (pink dashed curve). 



It should be noted that all measurements are in the 
free carrier concentration range set by parameters #2-1 
and 2-2. The 50 K data is very well fit with the same 
parameters as fit the 10 K data except 7 requires an in- 
crease from 7.5 to 13.5 cm^^. This suggests that only the 
free carrier damping rate changes over this temperature 
range. At higher temperatures, phonon parameters and 
effective free carrier concentration may also change. 

Figure [3] shows the transmission of circular polarized 
light and the complex Faraday angle (Faraday rotation 
and circular dichroism) together with a model calcula- 
tion with one set of parameters (#4). The green dash- 
dot curves were calculated using fit #1 parameters and 
clearly deviate strongly from the experimental results. 
The qualitative difference in shapes of the transmission 
curves for figures [2] and [31^ a) is also understood in terms 
of carrier concentration inhomogeneities. 

The sign of Wc in our magneto-optical experiments is 
determined by the sign of the carrier charge, the direction 
of magnetic field, the chirality of incident light, and other 
experimental conditions. To avoid any possible ambigu- 
ity, a boron doped Si wafer with an anti-reflection coat- 
ing in which the free carriers are holes was measured as 
a reference sample. The measured Faraday rotation for 
the Bi2Se3 flake has the same sign as for the Si:B wafer. 
However, the sign of Bi2Se3 carriers is opposite to holes 
because of a sign change of the Faraday response due to 
an etalon effect in the flake as we see from the model 
calculations. Therefore, the free carriers in our Bi2Se3 
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FIG. 4. (Color online) Reflectivity of a Bi2Se3 crystal with 
A^i, — 8.5 X 10^*. Inset: model calculations for the phonon 
line shape in magnetic field using data of Ref. @. Deviations 
in line shape and an apparent phonon frequency shift from 
the B=0 T (solid black) curve is due to the differences in 
the optical responses of left (dashed red) and right (dashed 
blue) circularly polarized light at fields of 8 T. Linearly or 
unpolarized light is an equal admixture of the two responses 
(solid green). 



FIG. 3. (Color online) (a) Transmission of circular polarized 
laser light at 42.3 cm~^, (b) real and (c) imaginary parts of 
tan^F {Of — complex Faraday angle). 



sample are electrons. 

These results also bear on electron interaction effects. 
While the Coulorab interaction is weak because of the 
large static screening (eo « 100), a significant electron- 
phonon interaction may be expected in this strongly po- 
lar material. In either case, the carriers scattering rates 
and effective masses reflect the significance of the many 
body interaction effects. Electron-electron scattering is 
expected to have T'^+uj'^ dependence from Coulomb scat- 
tering and an onset of scattering for w > wlo due to LO 
phonon emission. In addition, the optically measured 
cyclotron mass differs, in principle, from the value mea- 
sured from quantum oscillations measurements which in- 
cludes the self energy corrections to the band mass due 
to interactions. For the optical measurement of cyclotron 
resonance, the Kohn theorem insures that the resonance 
is controlled by the band mass provided that the electro- 
dynamics is sufficiently local. ^^'^^ In the Faraday geome- 
try of our experiments, the non locality leads to a Doppler 
shift of the cyclotron resonance at uj = ujc — qvp, where q 
is a characteristic electrodynamic wave vector.— For this 
system q « Wp/c « 10^° << Uc/vp, so that the expected 
shifts are very small. Therefore, the observed cyclotron 
mass should correspond to the band mass. The optically 
measured rricb = 0.15±0.01. The cyclotron mass deduced 
from SdH measurements on samples from the same batch 



was msdH = 0.15 ± .01. The near equality of these two 
values implies that mass enhancement from many body 
interactions is very small in this material. 

In figure 31 we show the room temperature refiectivity 
spectrum of a crystal with a higher free carrier concen- 
tration (iVb = 8.5 X 10^*) than the flake. The spectrum 
shows a plasma edge at 500 cm"-'^. The plasma edge fea- 
ture is very sensitive to the electron relaxation rate at 
this frequency. The electronic relaxation rate implied by 
this spectrum is 7b = 58 cm~^. It is noteworthy that 
this relaxation rate is relatively small considering that 
the plasma edge is significantly above the LO phonon fre- 
quency which is lulo ~ 150 cm~^. This relaxation rate 
is comparable to the relaxation rate measured at lower 
frequencies at room temperature so that LO phonon scat- 
tering does not add significantly to the relaxation rate. 

Therefore, from these observations, we see that the 
electron-phonon interaction effects are too small to be 
observed on the scale of the uncertainty of the effective 
mass and the impurity scattering rate. This surprising 
conclusion may result from the free carrier screening of 
the electron-phonon interaction. 

It was argued in Ref. 5 that the electron-phonon in- 
teraction is strong based in part on the anomalous line 
shape of the lowest infrared phonon. In fig. |4] (inset), we 
show results of our model calculation for this phonon in 
magnetic field at low temperatures. The line shape of the 
phonon in magnetic field is an average of two very dif- 
ferent reflectivities and R~ which results in apparent 
softening and broadening of the phonon line. Our trans- 
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mission data presented in fig. [T] also show insensitivity of 
the lowest phonon to applied magnetic field. Thus, the 
observations of Ref. @| on this phonon at least partially 
originate in cyclotron resonance physics rather than in 
electron-phonon and magneto-electric couplings. 

In conclusion, we have measured the bulk cyclotron 
mass for a low doped Bi2Se3 single crystal which is close 
to the mass observed in SdH measurements. Free carri- 
ers are electrons as we determined from the sign of the 
Faraday angle. We have shown that an effective Drude- 
Lorentz model well describes the magneto-optical data 
which is dominated by etalon effects in the slab sample. 



No signature of strong electron-phonon interaction is ob- 
served. 
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